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Abstract. We perform a multi-band statistical analysis of core-dominated superluminal active galactic nuclei (AGN) 
detected with Fermi Large Area Telescope (LAT) . The detection rate of 7-ray jets is found to be high for optically 
bright AGN. There is a significant correlation between the 7-ray luminosity and the optical nuclear and radio 
(15 GHz) luminosities of AGN. We report a well defined positive correlation between the 7-ray luminosity and the 
radio- loudness for quasars and BL Lacertae type objects (BL Lacs). The slope of the best-fit line is significantly 
different for quasars and BL Lacs. The relations between the optical and radio luminosities and the 7-ray loudness 
are also examined, showing a different behavior for the populations of quasars and BL Lacs. Statistical results 
suggest that the 7-ray, optical and radio emission is generated at different locations and velocity regimes along the 
parsec-scale jet. 



1. Introduction 

Extreme optical variability was a clear signature for some 
blazars detected at 7-ray energies during the EGRET era 
(Fitchel et al. 1994). Using the Fermi/LAT improved sen- 
sitivity, recent multiwavelength variability studies on indi- 
vidual sources have confirmed the tight connection among 
the 7-ray and optical variable emission (e.g. 3C 273 in 
Abdo et al. I2010ap . So far there was no attempt to com- 
pare optical continuum properties among 7-ray blazars. In 
this manuscript we take advantage of the unprecedented 
sensitivity provided by the first bright source catalogue 
(IFGL; Abdo et al. I2010bp to study the muhi-band cor- 
relations between optical, radio, and 7-ray emission in su- 
perluminal AGN detected by Fermi/LAT during its first 
11 months of operation. We adopt a flat cosmology with 
Ho = 71 km s-i Mpc-\ 17„ = 0.27, and = 0.73. 

2. The sample of AGN detected with Fermi/LAT 

The MOJAVE-2 sample consists of about 290 core- 
dominated AGN from the MOJAVE (Monitoring of Jets 
in AGN with VLBA Experiments) program (Lister et al. 
2009). The sample of AGN studied in this manuscript in- 
cludes 250 core-dominated AGN (Kovalev et al. 2005) 
from the MOJAVE-2 sample. Most of these AGN are 
currently monitored with the Very Large Baseline Array 
(VLBA) at 15 GHz. The IFGL catalogue includes about 
190 sources from the MOJAVE-2 sample. Optical nuclear 
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Fig. 1. Distributions of optical nuclear fluxes for AGN 
from the MOJAVE sample (full hue), 102 AGN de- 
tected as the gamma-ray source (filled histogram) and 
non-detected with Fermi/LAT during the IFGL period 
(shaded histogram). 



fluxes and redshift measurements were available for a sam- 
ple of 102 MOJAVE sources (Arshakian et al. I2010ap 
identified by Fermi/LAT (hereafter, M-IFGL sample). 
Seventy six out of 102 sources in the M-IFGL sample are 
part of the statistically complete MOJAVE- 1 sample. We 
will refer to these 76 sources as the Ml-IFGL sample. 
The Kolmogorov-Smirnov (K-S) test shows that the dis- 
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Fig. 2. Radio luminosity at 15 GHz against optical nuclear 
luminosity at 5100 A. 



tribution of redshifts in the MO J AVE- 1 sample and the 
Ml-IFGL sample are drawn from the same parent popu- 
lation indicating that the Ml-IFGL sample is not biased 
by redshift. 

The M-IFGL sample consists of 76 quasars, 24 
BL Lacs, and two radio galaxies which are excluded 
from further statistical tests. We use the non-parametric 
Kendall's r test to analyze correlations between indepen- 
dent variables and the partial Kendall's r test to account 
for the common dependence on redshift in the correla- 
tions between luminosities. Throughout the paper we as- 
sume a correlation to be significant if a chance probability 
P < 0.05. 

3. Interplay between 7-ray, optical, and radio 
properties of superluminal AGN 

The distribution of optical nuclear fluxes for the MO JAVE 
AGN detected and non-detected with Fermi /L AT is pre- 
sented in Figure 1. At first glance, the Fermi /L AT de- 
tected sources seem to have higher optical fiuxes than 
those with no-detection. This is further supported by the 
K-S statistical test: there is a significant difference (at 
a confidence level of 99.9 %) between the distributions 
of optical fiuxes of AGN detected and non-detected by 
Fermi /LAT. When using the Ml-IFGL sample, the dif- 
ference is significant at a confidence level of 96 %. We con- 
clude that the detection rate of 7-ray AGN is high for 
sources having high optical nuclear fluxes. 

We derive the optical nuclear luminosities (Lopt), to- 
tal VLB A luminosities (ivLBA), and the rest-frame radio- 
loudness {R oc ^vlba /<5'opt ) using the fluxes given in 
Arshakian et al. (|2010ap . 7-ray luminosities were com- 
puted using the Eq. (1) in Ghisellini et al. (2009), where 
S J {1^1,1^2) is the energy flux between 0.1 GeV and 100 GeV 
from the IFGL catalogue. Note that the 7-ray, optical, and 
radio luminosities are estimated from non-simultaneous 
observations. 
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Fig. 3. Optical nuclear luminosity against 7-ray luminos- 
ity (top panel), and radio (15 GHz) luminosity against 7- 
ray luminosity (bottom panel) for AGN from the M-IFGL 
sample. Labels in the top panel denote the corresponding 
population. 



Optical-Radio emission. Arshakian et al. (j2010ap 
found a positive correlation in the Lopt — Lvlba rela- 
tion plane for AGN from the MOJAVE-1 sample. They 
concluded that the correlation is due to the population 
of quasars and that the optical emission is non-thermal 
and generated in the parsec-scale jet. This is supported by 
studies of individual radio galaxies, 3C 390.3 and 3C120 
(Arshakian et al. I2010bl and Leon-Tavares et al. I2010p for 
which the link between optical continuum variability and 
kinematics of the parsec-scale jet was found. It was in- 
terpreted in terms of optical continuum flares generated 
at subparsec-scales in the innermost part of a relativis- 
tic jet rather than in the accretion disk. We confirm the 
Lopt — LvLBA positive correlation for a larger sample of 
M-IFGL quasars as well as no-correlation for BL Lacs (see 
Figure 2 and Tables 1 and 2). 
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simultaneous optical/7-ray observations and stronger vari- 
ability in the optical regime than that in the radio, and/or 
wider range of Dopplcr factors in the optical regime com- 
pared to the range of Doppler factors of the jet at 15 GHz, 
if the bulk of optical emission is generated in the relativis- 
tic jet and it is Doppler boosted. 

We report a significant positive correlation (> 99%) 
between and radio-loudness for quasars and BL Lacs 
(see Figure 4 and Table 1). The solid line in Figure 5 
represents the best fit to the data, 



logL^ = (0. 95 ±0.11) log i?+ (42. 76 ±0.25). 



(1) 



logR 

Fig. 4. 7-ray luminosity [L^) versus radio-loudness (i?) 
for M-IFGL AGN. The solid line represents the ordinary 
least-square fit to the data. 



The Lj (X R relation suggests that the strong 7-ray jets 
have progressively high Doppler factors (or faster speeds) 
in the radio domain compared to those in the optical 
regime. For quasars, the regression line is fitted by 



\ogL^ = (0.33 ±0.14) logi?+ (45.8 ±0.32), 
while for BL Lacs the best fit is, 



(2) 



Gamma-ray~Optical emission. We find a positive cor- 
relation between and Lopt (Figure 3). The correlation 
is significant for quasars from the M-IFGL and Ml-IFGL 
samples and it is stronger for Ml-IFGL quasars (Tables 1 
and 2), suggesting for a single production mechanism for 
7-ray and optical nuclear emission. 

Gamma-ray-Radio emission. Kovalev et al. (2009) re- 
ported a significant correlation between 7-ray and radio 
VLBA (8 GHz) emission for sample of -30 AGN. This cor- 
relation holds at high confidence level (> 99.9 %; Figure 3) 
for non-simultaneous measurements and a larger sample 
(M-IFGL) of AGN. This suggests that the powers aver- 
aged over the long time scales are correlated and, hence, 
the Doppler-factors of the parsec-scale jet in the gamma 
and radio domains are not changing substantially on a 
timescale of a few years. Pushkarev et al. (this proceed- 
ings) found that the 7-ray emission leads the radio emis- 
sion of the parsec-scale jet at 15 GHz with time delay of 
few months. They interpreted the observed time lag as 
a result of synchrotron opacity in the jet: the radio and 
7-ray emission are generated in the same region (pertur- 
bation in the jet?) and become observable with some time 
delay due to the opacity effects. If this scenario is correct 
then the variable optical emission is also generated in the 
perturbation moving upstream the jet, and we should ex- 
pect that the optical emission leads the radio emission and 
delays with respect to the 7-ray emission. 

The correlation in the £vlba — relation plane is 
apparently stronger than that in the Lopt — L-y (Figure 3). 
One may think that the corrected optical nuclear emis- 
sion of some AGN (Arshakian et al. I2010ap might be con- 
taminated by contribution of a stellar component thus 
causing the large dispersion in the Lopt — Lj diagram. 
Contamination should be stronger in radio galaxies and 
weaker in quasars and BL Lacs for which the contribution 
of optical nuclear emission is dominant. The large disper- 
sion in the Lopt — relation plane can be due to non- 



log L^ = (0.98 ±0.21) log i?+ (42.29 ±0.44). 



(3) 



It is noticeable that the slope derived for quasars is shal- 
lower than the slope fitted for BL Lacs. The ~ R cor- 
relation is still significant for all AGN (> 99%), quasars 
and BL Lacs (99%) of the Ml-IFGL sample (Table 2). 
The best-fit parameters for the later sample are almost 
unchanged. 

Table 1. Kendall's r correlation analysis between emis- 
sion characteristics of AGN from the M-IFGL sample. 
Al and A2 are the independent variables for which the 
Kendalls r correlation analysis is performed, r is the cor- 
relation coefficient, and P is the probability of a chance 
correlation. The correlations are considered to be signif- 
icant if the chance probability P < 0.05 (or confidence 
level > 95%). 



All 


Quasars BL Lacs 


Al A2 T P 


T P T P 


.^VLBI Lopt 0.2 7e-3 
L-, Lapt 0.1 4c-2 
L~, ivLBA 0.4 6e-8 
L-, R 0.4 3o-9 
Lapt 5-y/SvLBA -0.3 lc-5 


0.2 2o-2 0.1 0.4 
0.1 5g-2 0.1 0.3 
0.3 7c-5 0.5 7o-4 
0.2 8c-3 0.4 I0-4 
0.2 4o-2 0.5 I0-3 


Table 2. Kendall's r correlation analysis between emis- 
sion characteristics of AGN from the Ml-IFGL sample. 


All 


Quasars BL Lacs 


Al A2 T P 


T P T P 



Lvi^BA 
L-, 
L-, 
L-, 

Lapt 



Lapt 
Lapt 
i'VLBA 
R 

S-y/ SvLBA 



0.2 
0.2 
0.4 
0.4 
-0.3 



le-2 
2c-2 
lc-6 
3c-6 
lc-4 



0.2 

0.2 
0.3 
0.2 
0.2 



2e-2 

lG-2 

lc-4 
lc-2 
4c-2 



0.1 
0.1 
0.4 
0.6 
0.5 



0.6 
0.7 
5c-2 
lc-2 
2e-3 



Gamma-ray loudness. We define the 7-optical loud- 
ness and 7-radio loudness as the S^/S opt and S'^/S'vlba, 
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Fig. 5. Radio VLB A luminosity against 7-optical loudness (left panel), and optical luminosity against 7-radio loudness 
(right-panel) for the M-IFGL quasars and BL Lacs. The solid line represents the ordinary least-square fit to the data. 



respectively. We find that the radio luminosity is inde- 
pendent of 7-optical loudness (Figure 5, left panel). On 
the other hand, there is a negative correlation in the 
Lopt — Sj/Syi,ba relation plane (Figure 5, right panel), 
which is much stronger for BL Lacs (confidence level ^ 
99.9 %; Tables 1 and 2). For the M-IFGL AGN, the best-fit 
regression line is represented by Lopt oc (S'^/S'vlba)"^'^''- 
We interpret this relation as the optically weaker jets to 
have higher Dopplcr factors (or Lorentz factors if outflows 
radiating in optical and 7-ray bands have the same view- 
ing angle) in the gamma domain compared to those in the 
radio domain. The slope of the best-fit line is steeper for 
BL Lacs (-0.89±0.21) than for quasars (-0.24±0.1) and 
the diflFerence is statistically significant. We suggest that, 
on the average, the ratio between the Doppler factors in 
gamma and radio regimes varies slower for quasars than 
for BL Lacs, indicating for different velocity regimes in 
their jets. 

The significant differences found between quasars and 
BL Lacs are supported in recent studies (e.g. Ghisellini et 
al. 2009, Sambruna et al. 2010, Tornikoski et al., in this 
proceedings) suggesting the presence of different physical 
conditions along the jet in quasars and BL Lacs. This 
and other insights into the relationship between the 7- 
ray, radio, and optical emission will be pursued in a more 
detail in further studies. 

4. Summary 

Using the sample of 100 superluminal quasars and BL Lacs 
detected by Fermi/hAT, we investigate relations between 
their optical, radio, and 7-ray emission available from non- 
simultaneous observations. Our main results are summa- 
rized as follows: 

— The detection rate of 7-ray superluminal AGN is high 
for optically bright AGN. 



— The known positive correlation between Lopt and 
LvhBA holds for the M-IFGL quasars at a confi- 
dence level of 98%. The known correlation between 
7-ray flux and radio flux density (measured quasi- 
simultaneously) is also valid for non-simultaneous mea- 
surements of L^ and Lvlba for quasars. The — 
LvhBA correlation is significantly stronger than that 
in the L-y — Lopt relation plane. 

— There is a correlation between Lopt and L^ which ex- 
clusively holds for quasars. The correlation is signifi- 
cant for the M-IFGL quasars (at a confidence level of 
99 %) and marginally significant for quasars from the 
Ml-IFGL sample (c.l. 95%). 

— We report a statistically significant positive correla- 
tion (c.l. > 99 %) between 7-ray luminosity and radio- 
loudness for both, quasars and BL Lacs. The slope of 
the L-y — R relation is found to be steeper for the pop- 
ulation of BL Lacs. 

— We find that the radio luminosity at 15 GHz is in- 
dependent of the 7-radio loudness (Sj/Svlba) for 
quasars and BL Lacs. The 7-optical loudness (Sj/Sopt) 
and optical nuclear luminosity are negatively corre- 
lated for quasars (c.l. 96%) and BL Lacs (c.l. 99.9%) 
with the slope of the later being steeper. 
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